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Abstract 


In this paper we study the problem of the global existence (in time) of weak, entropic solutions 
to a system of three hyperbolic conservation laws, in one space dimension, for large initial data. The 
system models the dynamics of phase transitions in an isothermal fluid; in Lagrangian coordinates, the 
phase interfaces are represented as stationary contact discontinuities. We focus on the persistence of 
solutions consisting in three bulk phases separated by two interfaces. Under some stability conditions 
on the phase configuration and by a suitable front tracking algorithm we show that, if the BV-norm 
of the initial data is less than an explicit (large) threshold, then the Cauchy problem has global 
solutions. 
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1 Introduction 

This paper concludes a long analysis, begun in mm, concerning the global existence in time of weak 


entropy solutions for a system of conservation laws modeling phase transitions in a fluid. More precisely, 
the focus of the analysis is on the persistence of solutions with stationary interfaces. 


The system under consideration consists of three equations, namely, 



( 1 . 1 ) 


where the state variables {v,u,X) G P :=]0,-|-oo[xM x [0,1] denote the specific volume, the velocity and 
the mass-density fraction of the vapor in the fluid, respectively. System 0) is the conservative part of a 
more complex model first introduced in m, one of whose novelties lies in the pressure p, which depends 
not only on v but also on A. More precisely, the pressure is prescribed by the law 



( 1 . 2 ) 


where the function a is assumed to be strictly positive in [0,1]. For example, we can takep = {l+X)k/v 
for some positive constant k. 
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A related example of pressure law occurs in the theory of (non-isothermal) ionized gases. In that case, 
we have p = (1 + a)RT/mv, where i?, m and T are the universal gas constant, the molecular mass and 
the temperature, respectively, while a is the ionization degree; if a = 0 the gas is not ionized. In that 
model, however, a is usually assigned as a given function of both p and T by Saha’s law [T^] while, in 
the case under consideration, the mass-density fraction A is understood as an independent variable. 


System (I.ll is strictly hyperbolic in fl; the eigenvalues ±a/u are genuinely nonlinear while 0 is 
linearly degenerate. As a consequence, the field associated to this latter eigenvalue supports contact 
discontinuities, which are understood in the model as phase interfaces. 

The Cauchy problem for includes the initial data 


[voix),Uo{x), Xo{x)) , Voix) >v>0, X G 


(1.3) 


The global existence (in time) of solutions to the initial-value problem for any strictly hyperbolic system 
of conservation laws whose eigenvalues do not change type is well-known and can be proved either with the 
Glimm scheme or with a front tracking algorithm, see mm- The case of large initial data is a challenging 
problem and can be tackled only for special systems; then, the issue is to find classes of initial data with 
(large) total variation for which global solutions to the Cauchy problem exist. This topic was studied in 
m and m in the case of the isothermal, respectively isentropic, p-system, see also mm-. these results 
were extended in the seminal papers [MlIlT] to the case of nonisentropic gas dynamics. We also refer to 
m for the extension of Nishida’s result to the initial-value problem in Special Relativity. 

As far as (l.I)-(1.3) is concerned, a positive answer was first given in |3] and then in [S]. In particular, 
in the former paper an explicit threshold of the BV-norm of the initial data was provided in order to 
have the global existence of solutions. Moreover, both papers require a sort of balance of the BV-norm 
of the initial data: the larger the variation of (uo,Mo), the smaller the variation of Ao and vice-versa. It 
would be interesting to prove whether or not, for general BV-data, the solution exists globally in time. 
Motivated by the techniques introduced in [3], we considered in [T] the special case of initial data (1.31 
where 

if x < 0 , 

\r if X > 0 , 


Ao(x) = 


which models the dynamics of a two-phases fluid. In order to deal with this particular framework, an 
original Riemann solver was proposed by which we proved the global existence of solutions for a wide 
class of large initial data. Such results improved by far those of [3] when adapted to that setting. We 
briefly mention that the problem of smooth perturbations of Riemann-type solutions has been studied by 
many authors, see for instance |15) . 

The next step regarded the case of initial data with two phase interfaces, namely. 


{ Xi if X < a , 

Xm if a < X < 6 , (1.4) 

Xr if X > /> , 

where Xi,Xm,Xr are constant in [0,1] and a < 6 are the location of the waves; denote ag = a(A^), 
dm = a(Am), ttr = a{Xr). Clearly, this case is much more complicated than the previous one, because of 
the possible bouncing back and forward of the waves in the middle region [a, 6]. In [3] we answered in the 
positive to the above issue by assuming the condition am < minja^, a^}; if v is fixed, then the pressure in 
the middle region is lower than the pressure in the outside regions. We also recovered the results of [5 by 
passing to the limit when either Om or Om —t Or- To give a physical flavor to the problem, assume 

for a moment that the function a(A) is increasing with A, as is the interesting case in modeling. Then 
the condition above also means that the mass-density fraction of vapor inside [a, 6] is less than outside: 
for brevity, we address to this case as the drop case. We refer to lain [ 3 ] for more details on the model 
and references. 
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In this paper, we deal with the two remaining cases, namely, 

either Um > mui{ai,ar} or at < am < flrj 

since the case ai > am > <ir can be deduced by the latter. Reminding of the previous physical inter¬ 
pretation, we shall loosely address to these cases as the bubble case and the increasing-pressure case, 
respectively. As in [5], we introduce two special Riemann solvers. Both of them replace the standard 
non-physical waves 0 by waves defined through integral curves] this allows us to attach these waves to 
the phase waves. Then, following again the lines of [2], we introduce an “asymmetrical” Glimm functional 
F, in the sense that its interaction potential Q takes into account only certain shock waves approaching 
the phase waves while, on the contrary, all rarefaction waves approaching the phase waves are included. 
We point out that the definition of F differs in each of the three cases mentioned above. 

A key feature of both [2] and the current paper is the possible occurrence of a stability condition, 
depending on the case under consideration, to control interactions in the middle zone. More precisely, 
such a condition is needed both in the drop and in the increasing-pressure cases and imposes bounds to 
the strengths of the phase waves. However, in the bubble case we require no stability condition and, as 
a consequence, any phase wave is admitted. 

The plan of the paper is the following. The main result is stated in Sectionj^ In Section|^we introduce 
the Riemann solvers, the composite waves and, at last, we define the functional F] we also recall some 
background facts from Haig. Sections]^ andfocus on the bubble and on the monotone-pressure cases, 
respectively; in particular, we prove there that the functional F decreases at every interaction. The last 
Section deals with the two cases at the same time; by showing the convergence of the front tracking 
algorithm we conclude the proof of the main result. 


2 Main Results 


In this section we state the main results of this paper, which concerns the existence of solutions to the 
Cauchy problem ( |1.1| )-( [T3| ,( |1.4| ). Under the notation defined in the Introduction, we set 

^ ^ ^ - a, ^ ^ ^ ^ ar - am ^ ^2.1) 

Om F ag Or F Om 


The quantities rj and ( range over ] — 2, 2[; they are the strengths of the two contact discontinuities (see 
(3.4) below) that support the phase interfaces of the model, located at a; = a and x = 6, respectively. We 
denote by C,M.,TZ the three regions separated by 77 and C, in the (x,t)-plane, see Figure 

As already mentioned in the Introduction, we focus on the bubble case, which corresponds to consider 
77 > 0 and C < 0 : and on the increasing-pressure case, which corresponds to both 77 > 0 and C > 0 . 


V C 


C 


M 


n 


X — a X = 6 

Figure 1: The two phase waves with strengths 77 and C in the {x, t)-plane 
and the regions £,M,TZ. 


In order to state the existence theorem, we introduce some threshold functions. First, as in mm , we 
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define the strictly decreasing function 


:= 1 
1 + r V 


(^1 + vT+r^^ 


r e 


( 2 . 2 ) 


which plays a key role in the main results; we notice that limr_j,g+ /C(r) = +oo and limr_>.+oo = 0 - 
Moreover, as in [5], we need another function related to the stability of the two phase waves configuration, 
which differs from that in [2]. In the bubble case it is 


'Hb{\v\, ICI) 


4 

4- |?7CI 


max 



2+ICI 
2-ICI’ 


ICI 


2+hi ] 
2-NJ ’ 


(2.3) 


for (hi, hi) G Db := [0, 2[x [0, 2[ . As for the increasing-pressure case, the definition of the function He as 
well as that of its domain is more complicated though explicit; we refer to (5.23) and (5.22) below, 
respectively, and to Figure for a picture of Dc- We can immediately observe that in the bubble case 
the pair (h|, ICI) can vary inside the whole square [0, 2[x [0, 2[ , while in the increasing-pressure case they 
can cover only a portion of it. 


ICh 



Figure 2: The domain Dc in the (h|, |CI)-plane. 


We denote Po{x) = p [vo{x), Xo{x)). The following theorem states the global in time existence of 
solutions in the bubble and in the increasing-pressure case. Notice that the statement is the same in both 
cases. 


Theorem 2.1. Assume (1.2) and consider initial data (1.3),(1.4). Let D = Db and H = Hb in the 
bubble case, D = Dc and H = He in the increasing-pressure case. Moreover, assume that for rj, ^ as in 
(13, the pair (|? 7 |, |(C|) belongs to D. If 

1 


TV (log(po)) 


m\n{at,am,ar} 


TY[uc)<lC{HMX\)) 


(2.4) 


holds, then the Cauchy problem (1.1)-(1.3),(1.4) has a weak entropic solution {v,u,X) defined for t G 
[0, - 1 - 00 ). If rj = C = 0 the same conclusion holds with IC (H{\r]\,\f\)) replaced by -|-oo in \2.A\ . 

Moreover, the solution is valued in a compact set and {v{t, ■),u(t, ■)) G L°°([0, oo[; BV(IR)). 


The sub-level sets Sh = {(|f?|, |CI) € D : H{\ri\, |C|) < h}, h > 0, of the function H play an important 
role in condition (2.4), see for instance Figure in the bubble case. Indeed, condition (2.4) holds for 
every (I 77 I, |C|) € Sh if 


TV (log(po)) 


1 


min{af,Om,ar} 


TV {uc)<IC{h), 


since K. is decreasing. When h = 2, we have JC{2) = 21og(2 -|- •\/3)/3 and the case is particularly 
significative; let us consider, for example, the bubble framework. As in the drop case [5], the domain ^2 
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includes the segments [0, 2[ on each axis, but the 2-level set of T-Lb is no more the graph of the function 
C{\v\) = 2(2 - h|)/(2 -f |ry|), see Figure|3] 



Figure 3: Sets of level c of the function Hb'- cases c = 1, 2,3; the thin 
line is the curve ( = Cd^D- 


3 Preliminaries and Functionals 

In this section we collect some preliminary results from Ellllll]: focusing on the front tracking algorithm 
used to construct the approximate solutions. Moreover, we introduce some functionals needed to estimate 
the total variation of such solutions. 


First, we recall some basic notions. System (1.1) is strictly hyperbolic with two genuinely nonlinear 
characteristic fields (of family 1 and 3) and a linearly degenerate one of family 2. For i = 1,3, the f-th 
right Lax curves through the point t) = (D, ft. A) £ 17 are 


V I— >■ u + 2a{X)h{ei), Aj , n > 0 , 

where Si denotes the strength of an z-wave. 


1 


= o log T , 


£3 = ^ log 


and h is the function defined by 


h{e) = 


if e > 0, 
sinh e if e < 0 . 


(3.1) 


(3.2) 


(3.3) 


Rarefaction waves have positive strength, while shock waves have negative strength. The wave curve for 
the second characteristic field through (7 G 17 is given by 


A 


V - —,u, A 


aHX) 


X G [0, 1], 


and the strength of a 2-wave is 


(5 = 2 


o(A) — a(A) 
n(A) (^(A) 


(3.4) 
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We solve the Riemann problems by means of some Pre-Riemann solvers, which are introduced in the 
following proposition. We use the symbols ‘L’ and to denote Lax and Integral curves, respectively. 
For i = 1, 3 and 6 i S {L, /}, we define the functions 


0 ^ = 



if = L, 
if 6 , = I. 


(3.5) 


Proposition 3.1 (Pre-Riemann solvers). Fix 6 i G {L,/}, for i = 1,3. There exists a map Re^e^ ■ 
n X fl ^ Kx] — 2, 2[xK such that for any two states U- = (u_, u_, A_), U+ = (u+, m+, A+) G we have 

Re,eAU-,U+) = {euS,es), (3.6) 

where £3 represent waves of family 1,2,3, respectively, satisfying the following relations: 

e3-ei = ^logf^V a_0i(£i)+a+03(£3) = ^^^^, S = 2^±^. (3.7) 

2 \p_ / 2 a+ -I- a_ 


We denoted a± = a(A±), p± = p(v±, A±). 


We refer to [H Proposition 3.1] for a proof of the previous result. In a few words. Proposition 3.1 
states the existence of four Pre-Riemann solvers Rll, Rn, Ru and Ril, that prescribe how to solve a 
Riemann problem with Fwaves taken along (Lax or integral) j-curves, for i = 1, 3. In particular, Rix is 
the solver of mm that employs Lax curves. 


We use a front tracking algorithm [7] to build up the approximate solutions to (1.1)-(1.3|,(1.4). The 
first step in the construction consists in taking a sequence (Uo,Uo)i/gN of piecewise-constant functions 


with a finite number of jumps, that approximate the initial data (1.4) in the sense of Section 4]. We 
choose two parameters a = > t), p = p,y > 0 and proceed as follows. At time t = 0 we apply the 

solver Rix at each jump of the approximated initial data; we split rarefactions into a finite number of 
rarefaction shocks, each of size < a, whose speed equals the characteristic speed at the right state (see [H 
Section 4] for more details). Then, an approximate solution {v'',u^,X){-,t) is defined until the first time 
two wave fronts interact and a new Riemann problem arises. In the case of interaction between two fronts 
of families 1 or 3, we again use Rll and adopt the following strategy to approximate outgoing rarefaction 
waves: they are prolonged as a single discontinuity if they already existed before the interaction, otherwise 
they are split into a fan of waves as before. On the other hand, when solving an interaction of a wave of 
family 1 or 3 with a 2-wave <5, we possibly make use of a method that attaches certain reflected waves to 
(5; the outcome is a wave of the same family of the incoming one and a stationary composite wave, which 
is defined below. 


Definition 3.2 (Composite wave [2]). Consider two states U- = (u_,m_,A_) and U+ = {v+,u+, X+) 
of fl, with A_ ^ A_|_. The composite wave 60 = (<5q,^, 5q) connecting t/_ to U+ is the stationary wave 
defined by Jq = Rii(U-, U+). We write |5o| = If^ol + l^ol- 


Notice that in Definition 3.2 the waves i5q, i5o are given zero speed and i5o reduces to a 2-wave as long as 
<^0 = ^0 = 0- Hence, instead of dealing with 2-waves, we are left with composite waves belonging to a 
fictitious 0 -family. 

When a wave front of family 1 or 3 with strength Si interacts with a composite wave i5o = ((5g, 6 , (5q) at 
t > 0, we exploit two different procedures to solve the emerging Riemann problem of states U-,U+. We 
indicate by Ii{e){U) the integral curves of family i = 1,3, parametrized by e and of origin U G Tl] then, 
we define U- = Ii(Sq)(U-) and U+ = I 3 {— 6 q){U+). Here follows a description of the Riemann solvers; 
for a proof of the following facts see [H Proposition 3.7]. 


1. Accurate Riemann solver. If J^i] > p, then the solution is formed by waves £i,£o,£ 3 , where 
(£i,( 5 ,e 3 ) = Rll(U-,U+) and £0 = i5o. 
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^0 

(a) ( 6 ) 


Figure 4: A composite wave Jq in the (x, t)-plane: in (a) it is drawn 
as three parallel close lines, while the auxiliary picture ( 6 ) is used to 
determine the states in the interactions. 


2. Simplified Riemann solver. If |(5i| < p, then we distinguish case i = 1 and i = 3: 

i) for 1 = 1, the solution is formed by waves ei,£o such that (ei,5,£ 3 ) = R]j{U-,U+) and 
^0 = (^oi ^0 + ^3); 

ii) for 1 = 3, the solution is formed by waves £o ,£3 such that (£i,(5,£ 3 ) = Ril{U-,U+) and 

^0 = (^0 + ^1; ^o)- 


We emphasize that Qi — h in all cases and the following relations are verified: 


£3 - £1 = 


if z = 1 , 
if z = 3, 

sgn£i = sgndj. 


a_0i(£i) + a+03(£3) = 


a+0i(5i) 

a-QsiSs) 


sgn Ej = 


sgnJ • sgndi if z = 1 , 

—sgn 6 ■ sgn Si if z = 3. 


if z = 1 , 
if z = 3, 


(3.8) 

(3.9) 


We recall the interaction estimates [2 Lemma 5.2 and 5.4]. For z = 1, 3, we denote by Ei the strength 
of the transmitted wave and by Ej, j = 1,3, j ^ z, the strength of the reflected one (even in case of 
interaction with a composite wave treated by the Simplified solver). 


Lemma 3.3 (Interaction estimates). For the interaction between two waves at time t > 0 we have the 
following; let i,j = 1,3, j 7 ^ z. 


1. Assume that an i-wave Si interacts with a composite wave Sq = (Jg, <5, <5g). If jdij > p, we have 

1 


|£i - (5i| = |£j| <-|(5j(5| and |£g - ^g| = 0, 


(3.10) 


while, zfl^il < p, it holds 


|£i-(5i| = |£o-5o| = kil < 


where 




if Si <0 and either (i = 1, S > 0) or (i = 3, <5 < 0^ 


-ki(5| otherwise, 


Co = Co{p) = 


sinh p 
P 


(3.11) 


(3.12) 


2. If two waves ai and j3j of different families interact with each other, then. 


h\ = \o‘^\, \£j\ = \ff\- 


(3.13) 
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3. Assume that two waves at and j3i of the same family interact. If both on and /3i are shocks, then 
the reflected wave ej is a rarefaction, while the transmitted wave Si is a shock and satisfies 


\ei\ > max{|ai|, |/3i|} . 


(3.14) 


If cXi and Pi have different signs, e.g. ai < Q < Pi, then the reflected wave is a shock, both the 
amounts of shocks and rarefactions of the i-th family decrease across the interaction and one has 


kil < c{ai) ■ min{|ai|, \p,\} , c{z) = 


cosh z — 1 
cosh z + 1 


(3.15) 


About the term in Co in (3.12), we notice that C'o(p) > 1 for any p > 0 and Co{p) —t I"*" for p -n- O'*". 


Now, we introduce some functionals needed to prove the boundedness of the total variation of the 
approximate solutions. Using indices i,m,r to refer to waves in L,M.,IZ, respectively, we define 

i = l,3, 5i>0 

where ^ > 1 is a parameter to be determined. We denote 


i = l,3, 5i<0 
5i€C,M,-R. 


L = L^ + L^ + U, LO = kol + ICol 


and 


1 


L = L^ + L^+U =Y, k.|+ E 1^*1+ E lk| = ^TV(logp(t,-))-|r?o|-|Co|. 


i = l,3 
5, e£. 


i = l,3 
5i^M. 


i = l,3 


In the following sections we also introduce the interaction potential 


Q = Q^ + Q"^ + Q 


(3.16) 


and specify that differ in the bubble case and in the increasing-pressure case. The resulting 

functional 

F = L + L° + Q (3.17) 

is equivalent to the total variation of the approximate solutions and has an asymmetrical character in a 
double sense: firstly, it depends on the phases and, secondly, shocks and rarefactions play a different role. 
As for the latter, not only shocks are weighted by f (a procedure of [3] also exploited in mm), but the 
contributions of certain shock waves from Q are dropped. This is due to the fact that in the interaction of 


some shocks with a 0-wave we have that AL is already nonpositive; more precisely, from (3.8) it follows 


that such unnecessary waves are either 3-shocks interacting with a 2-wave (5 > 0 or 1-shocks interacting 
with a 2-wave (5 < 0. On the contrary, rarefaction waves are always counted in Q if they approach a 
phase wave. As shown in Figure the total variation of the solutions increases going towards the more 
liquid regions. 

In the potentials we insert some positive weights that keep track of the regions of 

provenience {C,A4,Tl) of the approaching waves and of the 0-wave approached {po or Co), see Figure]^ 
In the next sections we show that the functional F decreases across any interaction, under suitable 
conditions on the parameters C) p, which differ in the bubble case and in the increasing- 

pressure case. In the proof we fix mo > 0 and assume that, for i = 1,3, the size Si of any Ashock 
satisfies 

\5i\<mo. (3.18) 


Such parameter is necessary in our proof to relate the total variation of the initial data Uo,Po to the sizes 
of the phase waves. 






Figure 5: How the total variation varies for interactions with the phase 
waves in the bubble case (a) and in the increasing-pressure case (b). 



4 The bubble case 


In this section, we prove the decreasing of the functional F in 
phase-dependent interaction potentials are defined by 

= (if>| + iF^ICl) 5] Idal -f , 

(Jg > 0 Sg <0 

= l^aCI+C E l^3CI 

' 5]^>0 5j^<0 ' 53>0 53<0 

Si£M 5i<^M 5^eM 5^eM 

= ( k :^ h + ki \ c \) e E 

ai>o Si<o 

ilE-R Si&'R. 


(3.17) in the bubble case, where the 


As previously mentioned, the interaction potential Q in ( 3.16[ ) lacks certain shock waves: precisely, 3- 
shocks interacting with 770 and 1-shocks interacting with ^ 0 i see Figure In the next two propositions 
we list the conditions on the parameters and p needed for the decrease of F. 
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Proposition 4.1. Assume that at time t > 0 a wave Si, i = 1,3, interacts with one of the composite 
waves r]Q or (q. Then, Ai.F{t) < 0 provided that 


f>i, Kf;<K;, k^<kI 

Co{p) < 



Figure 8: Interactions of 1- and 3-waves with Co solved by means of the 
Accurate solver. The fronts carrying the composite waves are repre¬ 
sented by a single line. 


(4.1) 

(4.2) 

(4.3) 


Proof. Since the two cases give symmetric conditions, we only analyze interactions involving Co; see 
Figure]^ By (3.8 1 and (3.91 we have 


£3 — El — —(5i, 

£3 ~ £1 = <^3, 


|£i| - |iJi| — -|£3| , 
ksl - 11^31 = |£l| , 


if i = 1 , 
if i = 3 . 


i = 1 . 


If the interacting wave is a rarefaction, then by (3.10) we have 


^-1, 


AL + AL° = 


C|£3| + |£i| - l^il = (C- 1)|£3| < l^iCI 


if 1^11 > P, 

ksl + kil ~ l^il = 0 if l^il < Pi 

and AQ = Kif^\eir]\ — Kf\Sir]\ — Ar^|(5iC|- Therefore, since |ei| < |(5i| by the interaction estimates. 


AF < 


-K 


) I^ipI 


-b 


-Kl 




(Kf;-K:)\5,p\-Kl\5,C.\ 


if l^il > P, 
if l^il < P, 


which is nonpositive by (4.1 ( 3 ^ 4 . Instead, if the interacting wave is a shock, then in both the accurate and 
simplified cases we have AL + AL'^ = |£ 3 | -bCkil = ~(?~ 1 )|£ 3 | and AQ = f{K!^\ei\ — Kf^\Si\ 

As a consequence, AF < —(C — 1 )|£ 3 | + — Kf)\5i'q\, which is < 0 by ( |4.1[ )iq. 


i = 3. 


If the interacting wave is a rarefaction, then by the interaction estimates (3.10) AL + ALP = 


ksl + kil - = 2|ei| < I 53 CI and 


AQ = 


Kf;\e,p\-K^\5^C\ 
-^ri^sci 
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if l^sl > P, 
if l^sl < P- 



















Then, 


AF < 


1 + 


I<53CI 




l^sCI 


if li^sl > P, 

if I 1 J 3 I < P, 


which is nonpositive by (4.1)2, (4.2 )i. On the other hand, if the interacting wave is a shock, then 


ekil+C|e3|-?l^3l=2eki|<e|<53CI 


AL + AL° = 




kil + Cksl -= (^ + i)ki| < (? + i)-^I<53CI 


and 


Therefore, 


AQ = 




AF < 




l^sCI 

I<53CI 


if 1 ^ 3 ! > P, 
if |<53| < P- 

if l<^3| > P, 
if l^3| < P, 


which is nonpositive by (4.2)i, (4.3)i. 


if I<53| > P, 
if I<53| < P, 


□ 


Proposition 4.2. Consider the interaction at time t > 0 of two waves of the same family 1 or 3 and 

^(t) < 0 p 

1<?< 


assume (3.18). Then, AF(t) < 0 provided that 

1 


c(TOo) ' 


< 


' ^ if™ < '^ ^ 

ICI ’ - \v\ ’ 

(4.4) 

<|p| + ^c'ICI<?-i- 

(4.5) 


Proof. First, we consider the interactions taking piace in J\4, see Figure For brevity, we only deal 
with the case of interactions between two 3-waves 03 and /I 3 giving rise to £i and £3 (the 1-waves case is 
analogous). 



Figure 9: Interactions of 3-waves in M.. 


If both 03 and P 3 are shocks, then £i is a rarefaction by Lemma 3.3 and, as in [T] Proposition 5.8], we 
have 

AL-f |£i|(e-l) = 0, (4.6) 

for any ^ > 1. Moreover, we have 


AQ = Kf;\e^p\ - fK^\e,C\ < , 


AF < 


-(e-i) + i^rM kil 
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and F is non-increasing by ( 4 . 4 )i_ 3 . On the other hand, when the two interacting waves are of different 
type, for example 03 < 0 < /? 3 , as in [U Proposition 5.8], by (4.4)i one can deduce that 

AL + a^-l)\ei\<0. (4.7) 


3.3 


in both cases it holds 


If £3 is a rarefaction, then AQ = ^K^\eir]\ + (|£ 3 | — ^|a 3 | — |, 03 |) |C|; on the other hand, if £3 is a 

shock, AQ = ^K^\eir]\ + (^jesj — — l/^s]) jCI- Therefore, by Lemma 

AQ<CX™|£ir?|, + 


£1 


and F decreases by (4. 4 ) 1 ^ 3 . The analysis of the interactions between 1-waves requires symmetrically the 
condition < (C — 1)/ICI- 

Next, we analyze the case of interactions taking place in TZ, the case of interactions in C being 
analogous. By (4.4)i it is easy to verify that F decreases when two 1-waves interact. 



Figure 10: Interactions of 3-waves in TZ. 


Now, we consider the interactions between 3-waves (see Figure [^. When the interacting waves 
are both shocks we have (4.6 1 , while in the other two cases estimate (4.7) still holds under condition 
(4.4)i. If 03,/33 < 0 we have 

AQ = iF;|£i,7|+if^n£iCI, Ai^ = 


-(e-i) + if;H + ^cici kih 


while if, for example, 03 < 0 < /33 we have 


AF<i -{^-1 ) + k;\i^\ 


£il • 


Consequently, F is non-increasing by (4.5li. The condition K^\'q\ -|- K^\C,\ < C ~ 1 is required for the 
interactions occurring in the region C. □ 

Now, we can determine the order of choice of the parameters. To simplify the analysis, we can let 
and iF™ = since the final result does not change otherwise. Once 77 ,^ have been fixed, 
we choose in turn: mo, and AT™, and finally, we choose p (i.e. Co). First, notice that the 

conditions in (4.2) identify the set in the (AT™, iF™)-plane represented in Figure 11 Hence, by (4.2) we 
deduce 


AT™ > 


> l-b ^(l-biF^y) and 


AT™ > 


that imply 


T^m > 1 + ICI/2 
" - l-I^CI/4 


and 


i + if-M>i + M(i + i^-M), 

1 + I 77 I /2 


> 

< - I-I77CI/4 


(4.8) 


In particular, by replacing the inequality sign by equality in (4.8) we get the coordinates of the intersection 
point V between the two lines of Figure 11 Notice also that (4.8) implies ( 4 . 1 ) 2 . 
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Figure 11: Graphical representation of conditions (4.2) for I 77 I = 1/2 
and 1^1 = 3/2. 


Since we have chosen if™ = and = K^, conditions (4.5) imply ( 4 . 4 ) 2 , 3 . By ( 4 . 1)3 and (4.5)i, 
we get K^\t]\ + (^ — l)|Cl/2 < C ~ 1; which is equivalent to 

1^1 


K 


similarly, by ( 4 . 1)3 and ( 4 . 5)2 we get 


k: 


" 1 -ICI/2 
ICI 

1 - 1 ^ 1/2 


c 


<C-i; 


<C-i- 


(4.9) 


(4.10) 


By (4.8), (4.9) and (4.10) it follows that ^ must satisfy the inequality 

1 +ICI/2 hi l + hl/2 Id 


C > 1 + max 


1-ICI/2 l-hCI/4’ l-hl/2 l-hCI/4 


This condition must match with (4.4)i; then, recalling (2.3) we must require 

1 


1 + Hb(h|,|c|)<c< 


c(TOo) ’ 


(4.11) 


which is a condition that relates TOq to hi, |C|- When one of the phase waves tends to zero, 'Hf,(h|j ICI) 
tends to the other one and we completely recover the results of [T]. 

Summarizing, we choose the parameters as follows and keep strict inequalities for later need; let 
(hi, ICI) e Db be given. 


• First, we fix mo such that 


i + H,(hl,ICI) < 


1 


c{mo) 


(4.12) 


and take C in the interior of the interval given by (4.11). 


• In the (/C™,^P)"Plane we choose a point in the affine cone defined by (4.2) and sufficiently close 
to V ; moreover, we require 

1^1 ^ 1 T^m ICI 


K': 


<c-i 


and K, 


^ l-hl/2 


<C-1. 


1-ICI/2 

We choose and, then, by (4.131 we choose and such that 


?_^ ^ T^r ^ ^^ _ T^m \v\ Q-nrl 

2 ICI ICI ^ 


<K^n< ■ 

2 ^ hi ^ hi 


(4.13) 


(4.14) 
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Finally, we choose p such that (4.3) holds. 


Now, we can prove the global in time decreasing of the functional F. 


Proposition 4.3. Let mo > 0 satisfy ( |4.12[ ). Moreover, assume thatf^, ’ and p satisfy ( |4.11[ )-( |4.14 | 
and (4.3). Then, the following two statements hold. 


i) Local Decreasing. For any interaction at time t > 0 between two waves satisfying (3.18), it holds 

AF{t) <0. 


ii) Global Decreasing. If 

Z(0) < moc{mo) (4-15) 

and the approximate solution is defined in [0,r], then F{0) < mo, AFlt) < 0 for every t G (0,T] 


and (3.18) is satisfied. 


Proof. The first statement has been proved above. As for the second assertion, let us denote by Lf^’^ 
and the partial sums in due to z-rarefaction waves {iR) and z-shock waves (iS), respectively. 


By (4.4)2 ,3 we have 

F^{0) = L^{0) + Q™(0) < L™(0) (l + max{iF“|z7|, iF^"*|C|}) < . 


Moreover, from (4.5) it follows 


Then, 


F\0) < LUo) + LisiO) + (l + + i^s(O) (l + K^^\C\) < eL\0) , 

< TSiz(O) + LlsiO) + L[j,{0) (l + K;\r,\ + iF^HCl) + L[s{0) (l + K;\r,\) < . 

F(0) = F^iO) + F™(0) + F^O) < C^L(O). 


For a fixed t <T, suppose by induction that F{t) < mo and AF(t) < 0 for every 0 < t < t, interaction 
time. Then, the inequality AF{f) < 0 implies that 

F{t) < F(0) < e^L(O). 


Hence, by (4.15) the size of a shock 6 i (z = 1, 3) at time t satisfies 


1^*1 < < fL{0) < / x A(0) < mo 

f c(TOo) 


and (3.18) is verified. 


□ 


5 The increasing-pressure case 


In this section, we prove the decreasing of the functional F in (3.17) in the increasing-pressure case. As 
before, we first introduce the following interaction potentials 

Q^= (iF^M+AT^ICl) 5] |<53|, 

(53>0 

Ssez: 

= i^ii+^ E i<5ii)+^rici E 


Sg € A4 


5i>0 

Si€ 7 Z 


Si <0 

5i€TZ 
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The interaction potential Q in (3.16) lacks the 3-shocks interacting with rjo and with (q, see Figure 12 


The next proposition, which is analogous to Proposition 4.1 gives a first list of conditions that guarantee 
the decrease of F. 



Proposition 5.1. Assume that at time t > 0 a wave Si, i = 1,3, interaets with one of the composite 
waves 770 or Co- Then, i^Fit) < 0 provided that 


C>1, i^->i, 

(^ - <) H + (^r - ^c)ICI < 0, (1 - Kl)\C\ + + y) - 1^1 < 0, 

((e + i)y - ICI + C (^r(i + y ICl) - ^;) 1^1 < 0, (e + i)y - < o . 


Proof. Both in the case of interaction with 770 and Co, by (3.8), (3.9) we have 


£ 3 -ei = -<5i, |ei| - |(5i| = leal, if i = 1, 

£3 — ei = ^3, leal — l^al = —|£i|, ifi = 3. 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


However, in this case we have to treat separately the interactions with the phase waves, since the inter¬ 
action potential Q is not symmetric with respect to 770 and Co- 

Assume z = 1. If (5i is a rarefaction, then AL + AL° = 2 |£ 3 | < \Sir]\ and 


Interactions with 770 - 


AQ 


jK^\e3C\-Kff\S,p\ if|(5i|>p, 
\-Kff\d^v\ if|(5i|<p. 


Hence, 


AF < 


Then, AF < 0 by (5.3)2 and (5.1)2- 




\^lV\ 


1 - Kf 


\^iV\ 


On the other hand, if (5i is 


if l^il > P, 

if |5i| < p. 
a shock then 


AL + AL° 


2C|£3| < CI^ipI 

if l<5i| > P, 

(e + l)|£3|<y(C + l)IM 

if |<5i| < p 
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and AQ = —^K'^\5iiri\. Hence, 


AF < 


[ O - k\ l^i^l 


if l^ll > P: 

\Siv\, if|^i|<p, 


which is nonpositive by (5.1)2 and (5.4)2. 


Now, let 1 = 3. If S 3 is a rarefaction, then 




AL + AL°=-^'’’ ’ 2 


(C - l)|ei| < if|iJ3|>P 

0 if Idal < p 

and AQ = K^\s 3 (\ — K^\ 63 r]\ — K^\ 53 (^\. Then, since [eal < l^a] we have 

.e-1 


AF < 


- <) H + {K^ - ^c)IClJ I<53| if |(53| > P, 

{K^-K^^)\63C\-K\S3V\ if|<53|<p, 


and F decreases by (5.2)i. If Ja is a shock, then in any case we get AF = —(^— I)|ei| < 0 by (5.1)i. 
Interactions with Co- Assume i = 1. If (5i is a rarefaction, then AL + AL° = 2|ea| < |(5iC| and AQ 


K^\eir]\ — itr^|5i?7| ~ in both the Accurate and the Simplified case. Hence, by (5.2)2 we have 

(i-^c)ici + (^™(i + y)-^;) i^i 


AF < 

On the other hand, if is a shock, then 

^2C|e3|<CI<5iCI 


|<5i| <0. 


if l<5i| > p, 


AL + ALP = 


a 


(C + i)NI< + if|(5i|<P, 


and AQ = ^K'^\e 3 p\ - ^K';^\Sip\ - Thus, 


AF < 


(C +1) Y - ICI + c (^r(i + yI^O ■ 1^1 


if 1^11 > P, 


|(5i| if|5i|<p, 


which is nonpositive by (5.2)2 and (5.4)i. 


Assume i = 3. If (5a is a rarefaction, then 


and 


0 ,(e-l)ki|<^|^ 3 CI if|<5a|>P, 

AL + AL° = ^ 2 

0 if |(5a| < p, 

CF-|eiry|-F-|(5aC| if |(5a| > P, 


AQ = 


-KlSsCl 


if |<5a| < P- 
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Hence, 


AF < 


^ ^ I 7 ^™ C I ^ I f”! 


I<53CI if I<53| > p, 
if I<53| < P, 


which is nonpositive by (5.3)i. If, instead, ds is a shock, then AL + AL° = — (^ — l)|£i| and 

if l^sl > P, 


AQ = 


if |<53| < P- 


Then, 


AF < 


l)kil+^rki»7l if|^3l>P, 

-(C-l)|ei| if|^ 3 |<P, 


which is nonpositive by (5.1)i_3. □ 

As for an interaction between two waves of the same famiiy, Proposition|4.2|stiti hoids with the current 


functionai F. Therefore, the conditions required on the various parameters are ( |4.4[ ), (4.5). We omit the 
proof, since it can be carried out as above. 


Here, we make some comments on the conditions (5.1|-(5.4) and (4.4), (4.5); finaiiy, we estabiish the 
order in which we can choose the parameters. 

First, notice that (5.1)2 is implied by (5.3)2. Secondly, we can rewrite (5.2) as 

^bl + i^riCI<<H + ^clCh iF™(i + ^)H + |C|<if;H + iF[|CI- (5.5) 


Putting together (5.5)2 with (4.5)2 and (5.5li with (4.5)i, we have 

ICI 


a:: 


(i + y)h + ICI < + ifcnci < 1, 

^hi Ti^rici < ^>1 Ti^c'ici < c - 1 > 


then (4.4)2 ,3 are implied by (5.2) and (4.5). Moreover, by (5.3) we have 


ICI 


+ Y- 




that give the following lower bounds on and AT™: 

i + (C-i)ICI/4 


AT"* > 


i-CI^CI/4 


T^m > (C ~ 1) + Chi 
- 2(1-M/4)' 


Remark that (5.8) represents an affine cone in the (AT™, Ar™)-plane under the condition 

4 


c< 


\vC\ ' 


(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 


The vertex is the point whose coordinates are given by (5.9). Hence, A'™ and AT™ must be chosen in 
the non-empty intervals identified by (5.6), (5.9li and (5.7), (5.9)2, respectively. This means that (5.7), 
(5.9)2 give the condition 


C-l | I , iC^ll+Mi£i <y_i 

2 2(l-ChCI/4) ^ ’ 


(5.11) 
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while (5.6|, (5.9)i give the condition 


We introduce the notation \ri\ 
as, respectively, 


i-ei^ci/4 + + 


(5.12) 


= X, Id = y and ^ — l = z. Then, by (5.10) we rewrite (5.11) and (5.12) 


4 

(cy 

4 


(2 - x)z^ + y{x + 1) - (2 - a;) (^1 - ^) 


z + xy <Q, 


xy 


(4 - 2/) - 1 ^ + (l + 1)^ + y(l - ^) ^ 0- 


(5.13) 

(5.14) 


We also denote a{x,y) = xy{2 — x)/4, b{x,y) = y{x + l) — {2 — x) (l — xy/A), c{x,y) = xy, d{x,y) = a;y/4, 
e{x,y) = xy{A — y)/8 — 1 and f{x,y) = (l + y/2) x + y {l — xy/A) so that (5.13) and (5.14) become, 
respectively. 


a(a;, y)z‘^ + b{x, y)z + c{x, y) < 0 , 
d{x, y)z'^ + e{x, y)z + f{x, ?/) < 0 . 


(5.15) 

(5.16) 


Notice that the coefficients a, c, d, f are positive, e is negative and b may change sign. In order that each 
equations associated to (5.15) and (5.16) have distinct solutions, the discriminants 6^ — 4ac and — Adf 

(5.17) 


must be strictly positive. If 6 < 0, such solutions are positive. Thus, about (5.13) we require 

y{x + 1) — (2 — x) ^1 — + xyV‘2 — x < 0 , 


while about (5.14) we impose 

\xy 


{A-y)-l 


xy 


(i.|). + d-f) 


> 0 . 


(5.18) 


By a numerical comparison, we see that the set defined by (5.17) is included in that defined by (5.18). 
Under (5.17) and (5.18), we denote by 


zi, 2 {x,y) = 


(2 — x)(l — xy/A) — y{x + 1) ± y/ivix + 1) — (2 — x)(l — xy/A)]"^ — x^y‘^{2 — x) 


xy{2 — x)/2 

the two positive solutions of the equation associated to ( 5.15| ) and by 

1 - xy{A - y)/8 ± y/[xy{A - y)/8 - 1]^ - xy[{l + y/2)x + y(l - xy/A)] 


23,4(2:, 2 /) = 


xy/2 


the solutions of the equation associated to (5.16). Hence, by (5.10), (5.13) and (5.14) we get 

1 + ina,x{zi{x,y),Z 3 {x,y)} <^ < 1 + min j Z 2 {x,y), Zi{x,y), — - 1 

I xy 

Therefore, we can define the domain Dc represented in Figure as 

Dc = {(|??|, ICI) = {x,y) : ( |5.21[ ) holds} 

and the function 

■Hc(|?7l,ICI) = max{zi(|7?|,|C|),Z3(|77|,|C|)} . 


By (4.4)i we find the condition that relates nio to \r]\, |^|, i.e. 

1+77.(H,|C|)< ^ 


c(mo) ■ 


(5.19) 

(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 
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As a final remark, we notice that (5.4)i is equivalent to 


2? 




q^ici H<0. 


(5.25) 


Then, by taking p sufficiently small (since C'o(p) —^ 1 if p —>■ 0+) and C > !> (5.25) is implied by (5.2)2. 
For the choice of the parameters we proceed as follows. 

• We fix |? 7 |, |()| such that 

(5.26) 


(5.27) 


1 + max {zi{\r]l |CI), ^sdpl, ICI)} < 1 + min <{ zadpi, |C|), ^4(|??|, |CI), ^ f ’ 


Then, we fix ttIo such that (5.24) hold and, in turn, we choose ^ satisfying both (5.21) and 

i + ?7,(|7?|,|CI) <?< ^ 


c(TOo) ’ 


so that (4.4)i holds. 


• We choose AT™, AT™ such that (5.3) holds; in particular, we take (Tf™, AT™) sufficiently close to the 
vertex of the cone and satisfying (5.6)-(5.9), i.e. such that 


i + (e-i)|CI/4 „ e-i-ICI 

l-M/4 - ^ (1 + |CI/2)H’ 

Then, we choose A'^ = A'^, such that 


(^-l)|p|/2 + jF-|C| < ^ C-1 

H + ICI - " H + ICI’ 


(e-l)+el^l g-l-(e-l)|p|/2 

2(l-ChCI/4)- ^ ICI 


(5.28) 


iF;(,(l + |C|/2)|7;| + |C| 


H + ICI 




(5.29) 


Thus, (4.5) and (5.2) hold; hence, also ( 4 . 4 ) 2,3 are verified. 


• Finally, we choose p such that Co{p) verifies (5.4). 

In the next proposition, we claim the global in time decreasing of the functional F. The proof is 
omitted since it is analogous to that of Proposition |4.3| 


Proposition 5.2. Let mo > 0 satisfy (5.24). Moreover, assume thatf^, ’ and p satisfy (5.21)-(5.29) 
and (5.4). Then, the following two statements hold. 


• Local Decreasing. For any interaction at time t > 0 between two waves satisfying (3.18), it holds 

AF{t) <0. 

• Global Decreasing. If 

L{0) < moc{mo) (5.30) 

and the approximate solution is defined in [0,T], then F{0) < mo, AF{t) < 0 for every t G (0,T] 
and (3.18) is satisfied. 
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6 End of the Proof of Theorem 2.1 


In this last section, we conclude the proof of Theorem |2.1| and add some final comments. The proof fits 
in the general framework detailed in [2]; hence, we only outline the most important changes. 

First, as in Section 6] the front tracking algorithm used to construct the approximate solutions is 
well-defined and converges. Moreover, it is consistent in the sense that the total size of the non-physical 
waves carried by the composite waves vanishes with u. In brief, to estimate this quantity we use the 
notion of generation order, i.e. we attach an index fc > 1 to each wave generated in the construction. 
Then, according to each fc, we introduce suitable functionals Lfc, Qk, Ffe simply by referring the functionals 
L,Q,F to waves with order k. In particular, following the same steps as in mm, in both cases we can 
prove that 

j>k 

where ^ G ]0,1[ is either or ^c- More precisely, in the bubble case we find 


( 6 . 1 ) 


= max < 


l + iFpiCI l + 






1 + iF-hl l + iF^ICI 


2K^-1 


2K^ -1 ’ 1 + 2K^ ’ 1 -t 2Kl 




l + iF^H+iF^nCI 


Co 


i 

Co 


^ e - Co) ’ j ’ 

while in the increasing-pressure case we have 

+ + l + K^\ri\ l + if^dryl + ICI) l + K^^M + lCl) 


= max ■ 


2KJ^ - 1 


1 + 2K)! 


e 






?l^l/2 


ICI/2 


{K^ - l/2)|ry| + - iF-)|C| ’ (K^ - 1/2)|C| + [K^ - if™(l + |C|/2)]|r;| ’ 

_ ColCI/2 _ Co \ 

£,{Kl - a/2)|C| - K^{1 + Co|CI/2)]hl ’ - Co) j' 

In both cases, simple calculations show that /i < 1 because of our choice of keeping strict inequalities in 
the final parts of Sections and We exploit formula (6.11 to show that the total size of the composite 
waves tends to zero, as follows. We have 

[total size of composite waves] < 

< [size of composite waves of order > fc] -|- [size of composite waves of order < k] 


< ^ • F'i(O) -I- 7 ^ C'o(/o)(|?7| -|- Id) [number of fronts of order < fc] 


< ^ • Too 


P 


^ C'o(p)(|??| -I- Id) [number of fronts of order < fc] . 


( 6 . 2 ) 


Then, (6.2) is less than 1/:/ if we choose fc sufficiently large to have the first term less than l/{2v) and 
p = Pu{mo) small enough to have the second term less than l/(2z/). 


End of the proof of Theorem \2. 1\ As in the proof of [2 Theorem 2.1], in the bubble case (increasing- 
pressure case) by (2.4) and (4.15) ((5. 30]), respectively) we prove that 


F(0)< -TV(log(po)) 


1 


TV (uo) . 


Now, by (4.12) ((5.24), respectively) and 


2mm{ai,am,,ar} 
we look for an rUo satisfying both inequalities below: 
1 . 2 


(6.3) 


^(H,ld)< 


-1 = 


TV (log(po)) 


1 


min{ae,am,ar} 


c{mo) cosh Too — 1 

TV (mo) < 2moc{mo) =: z{mo) ■ 


=: w{mo ), (6.4) 

(6.5) 
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Recall from [2] that w{r) is strictly decreasing and z{r) is strictly increasing, for r G K"*". Moreover, we 
have JC{r) = z [w~^{r)), see (2.2). Hence, by (|2.4[) one can choose TOq such that (6.4), 


the cases. Therefore, we can conclude as in [7] and Theorem 2.1 is completely proved. 


(6.5) hold in both 

□ 


As in [2], we want to compare the results obtained here in the bubble casewith that of [3]. More 
precisely, we set a; = jTyj, j/ = |C| and we claim that 

'Hb{x,y)<x + y for 0 < a: + y < 1/2. (6.6) 

Since T-Lf, is a symmetric function of x and y, it suffices to verify that in the common domain it holds 

(2 + x)Ay 


(2 - x){4- xy) 


< x + y. 


(6.7) 


By simplifying expression (6.7), we find that it is equivalent to 

x'^y + xy^ — {2xy + 2y^ + 4a; + 8y) + 8 > 0 , 

which will be satisfied if a;y + + 2a: + 4y < 4. Since x < Xj^ — y^ this last inequality is verified if 

1 


^-yjy + y +2(^--yj+4y<4, 
that is when y < 6/5. Therefore, ( |6.6[ ) holds and, since K, is decreasing, we have 

lC{HbMAC\))>^{\ri\ + \C\) 


in the common domain |y| + |C| < 1/2. Hence, Theorem 2.1 improves [3l Theorem 2.2] in the bubble case. 
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